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Abstract 

We perform an analysis of general quark mass matrices in the general 
nearest neighbor interaction form. Excellent agreement with experiment 
is realized with this general texture, which is neither hermitian nor real- 
symmetric. We then propose a new class of quark mass matrices that 
contain no additional parameters other than the quark masses them- 
selves, and thus possess calculability, i.e., ensure a relationship between 
the six quark masses and four flavor-mixing parameters of the Standard 
Model. 
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I. INTRODUCTION 



A fundamental explanation of the flavor-mixing matrix, the fermion masses and their 
hierarchical structure persists to be one of the most challenging and outstanding problems 
of particle physics today. Within the standard model (SM) the fermion masses, the three 
flavor-mixing angles and the CP violating phase are free parameters and no relation exists 
among them. However, the expectation that "low-energy" quantities which can be computed 
in the SM should remain finite as the masses of intermediate particles go to infinity leads 
us to suspect that such a relation does hold. For example, AM(K° — Kf) diverges as 
m t — > oo. Since the contribution of the intermediate quark is always multiplied by a flavor- 
mixing matrix element, a relation between quark masses and flavor-mixing parameters could 
conspire to guarantee that the contribution to the low-energy quantity always remains finite. 

As an attempt to derive a relationship between the quark masses and flavor-mixing 
hierarchies, mass matrix ansatze based on flavor democracy with a suitable breaking so 
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as to allow mixing between the quarks of nearest kinship via nearest neighbor interactions 
(NNI) was suggested about two decades ago [1]. These early attempts are the first examples 
of "strict calculability" ; i.e., mass matrices such that all flavor-mixing parameters depend 
solely on, and are determined by, the quark masses. But the simple symmetric NNI texture 
leads to the experimentally violated inequality M top < 110 GeV, prompting consideration 
of a less restricted form for the mass matrices so as to still achieve calculability, yet be 
consistent with experiment [2]. 

It was later shown that the texture structure of these early ansatze for quark mass 
matrices, the texture structure of the NNI mass matrices, holds in general. Branco et al. 
have demonstrated that if one does not impose the assumption of hermiticity, then the NNI 
texture structure contains no physical assumptions [3]. For three fermion generations, one 
may consider without loss of generality quark mass matrices of the NNI form. This texture 
structure serves as the general starting point from which additional constraints on the mass 
matrix elements can be imposed in order to achieve calculability. 

In this paper we analyze general quark mass matrices in a modified NNI form to be 
described in the next section. This form is free from physical content, as it corresponds 
to a particular choice of weak basis in the right-handed chiral quark sector. We perform a 
numerical fit to the most recent measurement of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix, Vqkm- The results are in excellent agreement with the observables. Finally, we 
propose a new mass matrix ansatze based on our numerical results, to ensure some relation 
among the six quark mass masses and the four observable flavor-mixing parameters, thus 
reducing the number of free parameters in the SM. 

Our paper is organized as follows. In Section II we review the NNI form of mass matrices, 
pointing out a subtlety that has been overlooked in previous works. In Section III, we present 
the quark mass matrices to be analyzed as well as the resulting flavor-mixing matrix. We 
present the observables used in our fitting procedure as well as the resultant unitary triangle 
(UT) parameters. Guided by the numerics of the general case, we postulate quark mass 
matrices that ensure calculability. 

II. GENERAL NNI MASS MATRICES 

In a given generic gauge theory, flavor eigenstates do not necessarily coincide with the 
mass eigenstates which can be written in terms of combinations of flavor eigenstates. The 
mass matrix in a given electric charge charge sector is not necessarily diagonal in flavor 
space and involves the coupling between the left-handed and right-handed chiral states of 
different flavor. It can be shown that in the SM with three generations of fermions, one may 
perform a biunitary transformation of the quark mass matrices that leaves both the quark 
mass spectrum and the flavor-mixing parameters unchanged, where the new mass matrices 
in both the up-type and down-type quark sectors are of the NNI form and are in general 
neither hermitian nor symmetric [3]. We stress that this form is merely the consequence 
of choosing a particularly convenient weak basis that allows us to eliminate some of the 
redundant free parameters of the SM and has no physical consequences. 

Starting from initial completely general complex matrices M[ ud y a necessary condition 
for this biunitary transformation to yield the NNI form is the existence of a solution to the 
eigenvalue equation 
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(M^ + kM' d M'J).U t2 = XU, 2 



(1) 



where k is initially assumed to be an arbitrary complex constant. The NNI mass matrices 
are then given by 



M„ 



( A 12 e idl 

A 21 e id2 A 23 e W3 
V A 32 e ie * A 33 e i9 * 



(2) 



Ma 





B 21 e i9> * 




5126'' 


B 32 e l 



'i 

'< #33^' 



(3) 



M (M) = (4) 

where explicit forms for U and V( u ,d), corresponding to the left-handed and right-handed 
chiral quark rotations, respectively, are given in [3]. Note that U, the left-handed chiral 
quark rotation matrix, is common to both the up and down quark sectors. It is the freedom 
of performing an unitary transformation on the triplet of right-handed quark fields in both 
the up and down quark sectors, i.e. the freedom associated with the V( u ,d), without affecting 
the quark masses or the flavor-mixing parameters, that enables this texture structure to 
be possible while still remaining completely general. We can rewrite Eq. (1) for arbitrary 
complex k as 

All + A 2 23 + k(B 2 21 + B 2 23 ) = A 



Keai [cos(/3) sin (a) — sin(/9) cos(o;)] — k imag [sin(/3) sin(o;) + cos(/3) cos(a)] = 

where (3 = 9' 5 — 9' 3 and a = 6$ — 9 3 . The first equation exhibits the functional dependence 
of A on A; and provides no restriction on k. It may appear that when choosing a weak-basis 
that gives NNI forms for the quark mass matrices, one has two arbitrary degrees of freedom 
corresponding to k rea i and k imag . But in fact this two- fold degree of freedom does not exist; 
k must be either purely real or purely imaginary. 

For both k rea i and ki mag nonzero, we arrive at the inconsistency tan(/3) = — cot(/3). 
Inspection of Eq. (1) shows that if k is purely real, we are guaranteed a solution to (1) 
because the operator on the left hand side is Hermitian and has three linearly independent 
eigenvectors. However, this guarantee does not hold if k is purely imaginary. Therefore, 
we arrive at the condition k is purely real, which results in a — (3. Thus the number of 
parameters in our flavor- mixing matrix is reduced to eleven from twelve. 

The above mass matrices with their texture structure and the relationship a = (3 are the 
most general ones that exhibit the NNI texture and have the fewest number of independent 
parameters. At this stage no physical inputs have been introduced, and the resulting flavor- 
mixing matrix has five more parameters than are necessary to acheive strict calculability. 
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III. FLAVOR-MIXING MATRIX AND NUMERICAL FIT TO OBSERVABLES 



We use the rephasing freedom of the quark-fields to minimize the number of parameters 
entering into the quark flavor-mixing matrix. We can re-write Eq. (2) 

/ e i(*-94) o 0\ / A 12 \ / e *(fc-*3+fc) \ 
M u =\ e^-W 4 A a e^) = P (u)L M u P (u)R 

V 1 / V A 32 A 33 J \ J 

(5) 

with an analogous expression for M d , and where 0' 3 — 9' 5 = 3 — 0^. The flavor-mixing matrix is 
written in terms of the unitary matrices U( u )l and U^)l that diagonalize M U M\. and M d M d . 

U iu)L M u Mlu\ u)L = di&g(m 2 u ,m 2 c ,m 2 ) 

U (d)L M d M\u\ d)L = dmg(m 2 d ,m 2 ,m 2 b ) 
From the above expressions for M U:d , we see that 

M u Ml = P L M u M u T Pl (6) 

Because M U M U is a real symmetric matrix, it can be diagonalized by a real orthogonal 
matrix R u : 

R U M U MIR T U = diagK, ml m 2 ). (7) 
The flavor-mixing matrix can then be expressed as 

Vf-m = U( u ) L U {d)L = R u P {u)L P {d)L RT (8) 




Vf-m = h~ ^ ^ 1 \ \v d vZ v b \ (9) 



where the v U)C ,t are the normalized eigenvectors of M u mJ with eigenvalues m 2 u ,m 2 c and m 2 
and similarly for the v d , s ,b [4-6]. Note that only this single phase 9 = Q\ — 9^ — 9[ + 6' A enters 
the flavor-mixing matrix and is responsible for CP violation. Thus, it is sufficient for only 
one of the four non-vanishing elements coming from the second columns of M u and M d to 
be complex in order to provide a general parametrization of flavor-mixing and CP violation. 
For example, one may choose 6\ to be non-zero while all other phases vanish. At this stage, 
the mass matrices contain eleven independent parameters. 

We must specify the energy scale at which we are evaluating the mass matrices. We use 
the light quark masses [7,8] m u = 4.9 ±0.53 MeV, m d = 9.76 ±0.63 MeV, and m s = 187 ±16 
MeV, and the heavy quark masses m c = 1.467 ± 0.028 GeV, m b = 6.356 ± 0.08 GeV and 
m t = 339 ± 24 GeV all of which correspond to the masses at a modified minimal subtraction 
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(MS) renormalization point of 1 GeV. Choosing the energy scale to be 1 GeV we use 
the invariance of the trace, minor determinants and determinant of M U ^M^^ under a 

similarity transformation to relate the M u ^)ij to the mf in each quark sector. These six 
relations comprise the first six terms of our Xtot- 

In addition, we fit to the following V C km observables [9]: \V ud \ = 0.9740 ± 0.001, \V US \ = 
0.2196 ± 0.0023, \V cd \ = 0.224 ± 0.016, \V CS \ = 1.04 ± 0.16, = 0.0395 ± 0.0017, and 
= 0.08 ±0.02. 

Our method is to fit to the experimental data and compute the corresponding Xtot = 
Y^l=i Xi with one degree of freedom arising from twelve experimental values minus eleven 
independent input parameters. 

The following parameters yield a Xtot °f om y 0.180: 

Table 1 



Parameter 


Value [MeV] 


Parameter 


Value [Mev] 


A 12 


976.99 ± 0.9818 


B12 


44.149 ±0.1959 


A 21 


101.07 ± 0.975 


B21 


50.532 ± 1.406 


A 23 


-1465.5 ± 0.5671 


B23 


304.19 ±3.2851 


A 32 


-3.3811 xlO 5 ± 0.9675 


B32 


-3644.1 ±0.9721 


A33 


-24678 ± 0.9679 


B33 


-5203.4 ±0.9693 



The single phase 9 that minimizes Xtot * s -0-716 ± 0.0989 radians. The above parameters 
also make predictions concerning the mixing of the top quark as follows: |v fd |2 + |y* 6 |2 + |y tj[ |2 = 
0.9984, \V t * b V td \ = 0.009422 and j^j = 0.2479. The first two predictions above are in excellent 

agreement with the latest Particle Data Group values [9] , while our value of 1^4 is slightly 
higher than that predicted in [10,11]. 

The numerically determined mass matrices M u , M d from above display some interesting 
properties worth noting at this point. One can immediately see that they are far from being 
symmetric, as \An\ differs markedly in magnitude from \A 2 i\, etc. Of particular interest is 
the difference between the up and down quark sectors vis-a-vis the 3-3 element. Historically, 
mass matrix ansatze previously proposed in the spirit of calculability have written the 3-3 
element as a small pertubation about |m ti &|; i.e. as |m 3 — e|, where epsilon is some small 
parameter. Inspection of the above parameters shows that this description works fairly well 
for the down sector, but is not at all applicable to the up quark sector. In the up quark 
sector, it is the magnitude of the 3-2 element that differs only slightly from 1123 = nit- This 
numerical result and observed structure is a consequence of the top quark being so heavy. 

Using just the central values of the quark masses at the 1 GeV scale, the absolute values 
of the elements of the flavor-mixing matrix elements are: 



I VcKM I 



/ 0.974007 0.219610 .00316079 \ 

0.223856 0.975401 .0395023 
V .00942966 .0380247 0.999214 / 



To impose unitarity we need to go through one additional step after varying our input 
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parameters and fitting to the 
have in the matrices 



above experimental values. In our flavor-mixing matrix, we 



V 



V u 

v c 







T 


I) 


j and 


Vd 


Vs 


Vb 




V 1 


i 


i ) 



explicit expressions involving the mf, the six central values of the quark masses squared 
at the 1 GeV scale. The real symmetric matrices M U M U and M^M^ will not have these 
central values as their exact eigenvalues unless the fit is perfect; i.e. unless Y^=i xl — 0- 

To ensure unitarity, after having performed the fit to the experimental observations, we 
evaluate the eigenvalues of M u mJ and M d Mj to find the correct m 2 we should use in the vl 
and then re-evaluate our expression for V/_ m , thus ensuring an unitary flavor-mixing matrix. 

Although unitary, our derived flavor-mixing matrix contains unphysical phases that must 
be removed in order to perform an unitary triangle analysis. After finding the mass matrices 
parameters that yield the minimum x 2 , we put the flavor- mixing matrix into the standard 
CKM representation advocated by the Particle Data Group and then into an improved 
Wolfenstein parametrization from which we find the unitary triangle. 

The improved Wolfenstein parametrization can be obtained from the standard CKM 
representation with the following identifications [9,12]: 

A = S12 = sin# 12 , A\ 2 = s 23 = sin6» 23 , AX 3 (p - irj) = s 13 e~ i5 = sm6 13 e~ i5 



with modified Wolfenstein parameters 

p = p(l 



X -),rj = rj(l - ^-) 



2 n ' Iy 2 

With the above parameter values and using the true numerical eigenvalues of the resulting 
%<i)^>(i)! we arr i ve & t the following values for the flavor-mixing parameters: 

Table 2 



Parameter 


Value [rad] 


Parameter 


Value 


#13 


3.160065 x 10" 3 


V 


0.3444959 


#12 


0.22825323 


P 


0.00422775 


#23 


3.95070522 x 10~ 2 


a 


71.61961° 


5 


1.5585246 


P 


19.083502° 


A 


.22627623 


7 


89.296884° 


A 


.771433972 


J 


2.748573468 x 10~ 5 



J is the Jarlskog invariant [13,5], a rephasing invariant of the mixing matrix, and 
is given in the standard representation of the quark flavor-mixing matrix as J = 
Si2Si3S23Ci2C? 3 c 23 s ( 5 = 2.748573468 x lO" 5 . 

After imposing unitarity as described above, the absolute values of Vqkm become: 



CKM\ 



( 0.974058 0.226275 .00316006 \ 

0.226101 0.973303 .0394979 
V .00941607 .0384891 .999215 / 
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As alluded to in the introduction, one of the most appealing features of the NNI mass 
matrices is the large number of texture zeroes. One may take advantage of this texture 
structure to eliminate three of the five real parameters in each sector, expressing them 
solely in terms of the m 2 , up to a sign ambiguity associated with the square root branches. 
That is, we may use the invariance of the characteristic equation under a similiarity trans- 
formation to express ^23,^32, and A33 in terms of ml, ml, and mf, and similarly for 
the down quark sector. Choosing A± 2 and Ai 2 to be the as yet undetermined parame- 
ters in the up-quark mass matrix, and using the branch information revealed by Table 

1, one easily verifies the following relations: A 33 = - m ^X7i ' ^ 32 = -\/ b+y/ f^ , A 



1-23 



ml + m 2 c + ml - A\ 2 - A\ - A 2 2 - ^ggs where b = m 2 u + m 2 c +m 2 - 2A 2 2 - 

and c = A\ 2 {m\ + m\ + m 2 t - A\ 2 ) + m *™} m t - m 2 u m 2 c - m 2 u m 2 - m 2 c m 2 t . 

Similarly, for the down-quark sector, choosing B\ 2 and B 21 to be the parameters as 
yet undetermined by the invariance relations, one finds the following equalities: -B33 = 

B32 = _/SM, B 23 = ^m 2 +ml + m 2 - B\ 2 - B 2 21 - B 2 2 - where 

d = m 2 + m 2 s + m 2 - 2B 2 2 - ?§0< and e = B\ 2 (m 2 + m 2 + m 2 - B 2 2 ) + - m 2 m 2 s - 

11 11 

171 d m b — m s m b- 

Using the above mass matrices with their associated five total degrees of freedom corre- 
sponding to Ai 2 , A 2 i, B12, B 2 i and 9, and the six experimental values of IVc^mI observables 
as before, we obtain Table 3, with a xiot °f 2.182. This xiot * s no ^ nearly as good as in the 
previous case, but the simple explanation for this apparent decline in numerical confidence 
is just a reflection of the slight discrepancy between unitarity and the present experimental 
values. In the previous case, the flavor-mixing matrix could deviate from unitarity to satisfy 
the experimental constraints; i.e., \V CS \ having a central value that in itself violates unitarity 
did not do violence to the minimization of Xtot because the requirements of unitarity were 
not automatically satisfied. In the current situation, unitarity is automatically imposed, 
and so the increase in x 2 tot is to be expected. The predicted IVc^mI all lie within the ranges 
preferred by the Particle Data Group analysis, and as such, any set of parameters that 
fulfill this requirement of agreement with the ranges listed therein is deemed "acceptable" , 
regardless of xfot values. 

, .975244 .22111 .00314624 \ 
\V CK m\ = I -220939 .974488 .0394887 

.00924305 .0385204 .999215 / 

Implicit in this prejudicial "acceptance" is the belief that only three generations of quarks 
exist and that future experiments will bring the experimental values into closer agreement 
with the constraints imposed by unitarity. 
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Table 3 



Parameter 


Value 


A 12 


1229.3 ±32.38 


A 2 i 


139.82 ± .091 




45.872 ± 3.4573 


B 2 i 


48.797 ± .0433 


9 


-0.67806 ± .009657 



In addition, the standard representation parameters as well as the Wolfenstein parame- 
ters and Jarlskog invariant J are predicted to be: 



Table 4 



Parameter 


Value 


Parameter 


Value 


#13 


3.146249 x 10~ 3 


77 


.3515246 


#12 


0.2229535 


P 


.001394 


#23 


3.9499098 x 10~ 2 


a 


70.83438° 


5 


1.5668302 


(3 


19.39285° 


A 


0.2211108 


7 


89.77276° 


A 


.80770906 


J 


2.67698466 x 10~ 5 



Because the mass matrices squared have the correct eigenvalues, the free-parameters 
A12, A21, B12, B 2 i and # represent the general parametrisation at the fundamental mass ma- 
trices in general. To realize calculability, it is not enough to simply postulate some relation 
among this set of parameters independent of the quark masses, as such a relation would only 
serve to relate the flavor-mixing parameters among themselves, and not provide the sought 
after relation between quark masses and flavor-mixing parameters. 



IV. MAKING NEW MASS MATRIX ANSATZE 



The above mass matrices are completely general; to achieve calculability we must find 
expressions for A i2 , A 2 i, Bi 2 , B 2 i and 9 in terms of the rrii. Guided by the parameter values 
in Tables 1 and 3, we postulate the following mass matrices in terms of the mf. 



( ° 

A21 

V 



A 12 








\ 



-yjml + ml + m? - A\ 2 - A 2 21 - A\ 2 - 



A 21 



b+Vb 2 +4c 
2 



where A 12 = + yj*f) , A 2 i = ^m u m c and 

/ B 12 
M d = B 2 i 

v -v ^+vy+Ii 



A 12 A 2 







m\ + m 2 s + ml — B\ 2 



B 2 



21 



R2 
-°32 



a s b 
^12^21 



m b ^fm d m s 
B21 
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where Byi = ^m,dm s , B21 = + ^Jm^rrT s and b, c, d and e are defined as before. 

Lastly, the single phase 9 entering the flavor-mixing matrix is postulated to be 
We know that for only two quark generations there is no CP violating phase in the mixing 
matrix, so it is natural to expect, within the framework of calculability, that 9 will involve 
ratios of elements in the mass matrices that only involve the mixings of the third generation 
quarks. With these postulates for M u , M d and 9, V/_ m is found from Eqns. (6-9). This 
constitutes a new mass matrix ansatz that provides a calculable model of flavor-mixing and 
is in excellent agreement with the latest experimental results. 

This new ansatz predicts the following absolute values for (V^mI: 

/ .974427 .224682 .00307663 \ 
\V C km\ = .224513 .973672 .0394506 , 
V- 00923256 -0384783 .999217 / 

1^1 = 0.0779868 and j^j = 0.239942. The xtot fo und from the six experimental V C km 
measurements is found to be 5.2489. Because we have no free parameters, there are six 
degrees of freedom. Such a value for Xtot corresponds to ~ 70% confidence level. 

The standard representation parameters, Jarlskog invariant J, and the modified Wolfen- 
stein parameters are predicted to be: 



Table 5 



Parameter 


Value 


Parameter 


Value 


#13 


3.0766324 x 10~ 3 


A 


0.781473947 


#12 


0.226617779 


V 


0.3380146 


#23 


3.94622992 x 10~ 2 


P 


0.01469470 


S 


1.52735011 


a 


73.554460004° 


J 


2.652858 x 10~ 5 


P 


18.93482465° 


A 


0.2246832 


7 


87.51071528° 



The flavor-mixing observables predicted from this calculability ansatz in Table 5 are in 
excellent agreement with the general result of Table 2. 

The above mass matrices with their calculability property as well as the low Xtot are ver y 
compelling arguments in favor of calculability in the quark sector, but naturally one would 
like to uncover at least a glimpse of the more fundamental theory beyond the SM that is the 
source of this calculability. When one considers hermitian mass matrices, considerations of 
family permutation symmetry and its breaking in successive stages are suggestive explana- 
tions of the source of calculability, whereas with general NNI texture mass matrices, even 
though calculability still holds, it is not readily apparent what familiy symmetry, if any, 
is responsible. In conclusion, we have elucidated an important point in the construction 
of NNI weak-eigenstate quark basis that has eluded some previous authors. We have then 
performed an analysis of the experimental data to determine the mass matrix parameters, 
and have obtained values for the flavor-mixing and UT parameters that are in excellent 
agreement with previous analysis [12]. Finally, we have presented a new class of calculable 
mass matrices that also can explain all the experimental data and predict nearly identical 
values for the flavor-mixing and UT parameters as in the previous case. 
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